HEAT TRACE OF NON-LOCAL OPERATORS 
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ABSTRACT. This paper extends results of M. van den Berg on two-term asymptotics for the trace of 
Schodinger operators when the Laplacian is replaced by non-local (integral) operators corresponding to 
rotationally symmetric stable processes and other closely related Levy processes. 
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1. Introduction 

There is an extensive literature of trace formulae for heat kernels of the Laplacian and its Schrodinger 
perturbations in spectral and scattering theory. For instance, van den Berg [ ] and later Banuelos 
and Sa Barreto [5], explicitly compute several coefficients in the asymptotic expansion of the trace 
of the heat kernel of the Schrodinger operator —A + V as t J, for potentials V G 5(1^), the 
class of rapidly decaying functions at infinity. In particular, in [5] a general formula is obtained for 
these coefficients by using elementary Fourier transform methods. For applications of these results 
to problems in scattering theory, see [5] and references therein. The results in [ ] were extended by 
Donnelly [20] to certain compact Riemannian manifold. Heat asymptotic results have also been widely 
used in statistical mechanics, for details we refer the reader to the article [26] of E. H. Lieb on the 
second virial coefficient of a hard-sphere gas at low temperatures, and to the article [28] of M. D. 
Penrose, O. Penrose and G. Stell on the sticky spheres in quantum mechanics. We also refer the reader 
to Datchev and Hezari's overview article [ ] for various other related spectral asymptotic results and 
applications. 

For non-local (integral) operators which arise by replacing the Brownian motion with other Levy 
processes, many questions concerning spectral asymptotics are wide open at this point. These in- 
clude, for example, the Weyl two-term asymptotics for the spectral counting function of the fractional 
Laplacian with Dirichlet boundary conditions and the McKean-Singer [ 17] result involving the Euler 
characteristic of the domain in the third-term asymptotics of the trace of the heat semigroup. On the 
other hand, a two-term trace asymptotics involving the volume of the domain and the surface area of 
the boundary is proved in Banuelos and Kulczycki [3] (for C > domains) and in Banuelos, Kulczycki 
and Siudeja [4] (for Lipschitz domains) for the fractional Laplacian associated with symmetric sta- 
ble processes. These results are in parallel to the results of van den Berg [7] and Brown [10] for the 
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Laplacian. For further recent work on the Dirichlet case in domains of Euclidean space, see Frank and 
Geisinger, [22] and [23]. The purpose of this paper is to obtain analogues of the van den Berg result 
in [6] for the fractional Laplacian A"/ 2 associated with symmetric a-stable processes, < a < 2, 
in M. d and for other closely related non-local operators corresponding to sums of stable processes and 
the relativistic Brownian motion. These are all Levy processes which are obtained by subordination of 
Brownian motion. 



2. Stable processes, statement of results 

Before we state our results precisely, we recall the basic definitions and elementary notions related 
to stable processes. Let X t be the d-dimensional symmetric a-stable process of order a G (0,2] 

in M. d . The process X t has stationary independent increments and its transition density p[ a \x,y) = 

p[ a \x — y), t > 0, x, y G M. d , is determined by its Fourier transform (characteristic function) 



(2.1) = E(e* x <) = [ e**^ \y)dy, t > 0, £ € 

JR d 



i>d 



Denoting by P x and E x the probability and expectation, respectively, of this process starting at x we 
have that any Borel subset B C R d , x G R d , t > 0, 



P*(X t eB)= [ p[ a) (x-y)dy, 

JB 



where 

(2.2) pW( l)= * / e -*v-*rde. 

(2vr) a J Rd 
Here we can also write 

where r]^ 2 (s) is the density for the a/2-stable subordinator. While explicit formula for the transition 
density of symmetric a-stable processes are only available for a = 1 (the Cauchy process) and a = 2 
(the Brownian motion), these processes share many of the basic properties of Brownian motion. It also 
follows trivially from (2.3) that pf^ 1 (x) is radial, symmetric and deceasing in x. Thus exactly as in the 
Brownian motion case we have 

(2.4) p^ix) <t- d / a p[ a) (0). 

In fact, if we denote by the surface area of the unit sphere in W d we can compute p{^ (0) more 
explicitly. 

p(«)(0) = _L- I e -l«l> = He-'sG-Vds 
11 K ' (27T) d J Rd (27r) d a J 

{ ' (27r) d a ' 

Throughout this paper we will deal with radial transition functions and use the notation 

pf\ x ,y) =Pt \ x -v)- 
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In particular, p[ a \x, x) = p[ a \o). Of importance for us in this paper is the scaling property of these 
processes. More precisely, by a simple change of variables in (2.2), we see that these processes are 
self-similar with scaling 

(2.6) P ( t\x,y) = 

The transition densities p[ a ^ (x) satisfy the following well-known two sided inequality valid for all 

x G R d and t > 0, 



(2-7) (V rf /« A j-L^j < p^ (x) < C a , d (t- d /« A 



x\ d ^ 



where the constant C a ^ only depends on a and d. Here and throughout the paper we use the notation 
a f\b = min{a, 6} and a V b = max{a, 6} for any fl,iel. 

For rapidly decaying functions / G 5(IR d ), we have the semigroup of the stable processes defined 
via the the Fourier inversion formula by 



T t f(x) = E x [f(X t )]=E°[f(X t + x)] 

f(x + y)p t {dy) =p t * f{x) 



1 f e-<e-W a f(m- 

71") jR<i 



(27T) 

By differentiating this at t = we see that its infinitesimal generator is A Q / 2 in the sense that 
A a / 2 /(£) = —\C\ a f(0- This is a non-local operator such that for suitable test functions, including all 



functions in / G C^°(M ), we can define it as the principle value integral 



(2.8) A^f(x) = A d>a lim / f{X+ i V Lj {X) dy, 



where 



At - { 2 ' 



We will denote by Hq the fractional Laplacian operator A a > 2 , a G (0, 2] and by H its Schrodinger 
perturbation H = A Q / 2 + V, where V G L°°(]R d ). We let e~ tH and e~ tH ° be the associated heat 
semigroups and let pf and p[ a ^ denote their corresponding transition densities (heat kernels). In fact, 
Pt is just as in (2.1) and the Feynman-Kac formula gives 



(2-9) p?(x,y)=p\ a >(x,y)El y (e 
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where E x is the expectation with respect to the stable process (bridge) starting at x conditioned to be 



Tr(e- tH - e- tH ° 



(2.10) 



at y at time t. The main object of study in this paper is the trace difference 

- /tf <,.„)- 

= t- d l a p^ (0) I El x (e~ $ V ^ ds - l) dx, 

where p^\o) is the dimensional constant given by the right hand side of (2.5). 

Before we go further let us observe that this quantity is well denned for all t > 0, provided V G 
L°°(M. d ) n L 1 (R d ). Indeed, the elementary inequality \e z — 1| < \z\e^ z ^ immediately gives that 



EiJe-fonxs)ds_ 1 ) dx 



< e 



tl VI 



El 



\V(X a )\ds dx. 



However, 
(2.11) 



E\ 



^\V(X s )\ds) = £ EijV(X s )\ds 



p ( f\x,y)p ( fj s {y,x) 



\V(y)\dyds. 



/o Jm d pf\x,x) 
The Chapman-Kolmogorov equations and the fact that p[°^ (x, x) = p^ (0, 0) give that 

pi a \x,y)p[°^ s (y,x) 



-dx = 1 



and hence 
(2.12) 

It follows then that 
(2.13) 



' Rd Pt \ x i x ) 
j (j \V{X a )\d^dx = t\\V\\ 1 . 



|Tr(e 



o)| < t- d / Q+1 pi a) (0)||V-||ie 



illVI 



valid for all t > and all potentials V G L°°(M. d ) n L x ( 

The previous argument also shows that for all potentials V G L c 



(2.14) 



Tr( e -^-e-^)=pS Q) (0)^ 



("l) fc 



fc=i 



o 



t \ k 

V(X s )ds) dx, 



where the sum is absolutely convergent for all t > 0. 

We have the following two Theorems which parallel the results in van den Berg [ ] for the Laplacian. 
Note however, that here (ii) in Theorem 2.1 (as well as Theorem 2.2) is more general as we do not 
assume that the potential is nonnegative. 



Theorem 2.1. (i) Let V : 



(-oo,0], V G L°°(R d ) n L 1 (R d ). Then for all t> 



(2.15) J#(0)i||V||i < Tr(e- tH - e"^ ) < p[ a) (0) [ t\\V\\i + -^||V||i||V| 



„t||VI 
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In particular 



(2.16) 



Tr(e 



pi a \0) (t\\V\\ 1 + 0(t 2 )) 
t~ d / a p[ a \0) (t\\V\\ 1 + 0(t 2 )), 



as t 4- 0. 

(ii) If we only assume that V G L°°(R d ) n L 1 (R d ), tfien/or a// i > 0, 



(2.17) 



Tr(e-^ -e- ulo )+pj a) (0)i / V(x)dx < pf ; (0)C« 2 ||V r ||i||V r ||ooe t||v||o °, 
/or rome universal constant C. From this we conclude that 
(2.18) Tr(e" tH - e"^ ) = p[ a) (0) ^-t ^ ^ + 0(t 2 )^) , 

as t J, 0. 

Theorem 2.2. Suppose V G L°°(R d ) n L 1 (]R c( ) atic? if jj a/so uniformly Holder continuous of 
order 7 (there exists a constant M G (0, 00) smc/z f/iaf | V(a;) — V(y)| < M|x — y| 7 ,/or all 1,1/6 
wjY/i < 7 < a A 1, whenever < a < 1, and with < 7 < 1, whenever 1 < a < 2. Then for all 
t > 0, 



(Tr{t 



AH 



)) +p[ a) {0)t I V(x)dx-p\ a> {0)^t 2 I \V(x)\ 2 dx 



(«), 



1 



<c7 a , 7 , rf ||y|| lP ! a) (o)(|y||Le i ii y ii-t 3 + ^ Q+2 ), 



Tr(e' tH - e~ m °: 



pi a) {0)(-t [ V{x)dx + \t 2 [ \V{x)\ 2 dx + 0{t< /a+2 ) 



(2.19) 

where the constant C a ^ y d depends only on a, 7 and d. In particular, 
(2.20) 
as 1 10. 

Theorems 2.1 and 2.2 are proved in §2. In §3, we provide extensions to other non-local operators. 

3. Proof of Theorems 2.1 and 2.2 

Proof of Theorem 2.1. Setting 



/ V{X s )ds, and b = t\\V\\ 00 , 
Jo 



we observe that — b < a < 0. We then use the following elementary inequality from [6] 

1 



(3.1) 



( 1 + ^be b 



to conclude that 



J V{X s )ds < (e-tiv(.x.)d, - l) < -J V(X s )ds 



1 + 2 £ " V " l|o ° e 



film 



Taking expectations of both sides of this inequality with respect to E x x and then integrating on M. d 
with respect to x, it follows exactly as in the derivation of the general bound (2.13), that (2.15) holds. 
This concludes the proof of (i) in Theorem 2. 1 . 
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For (ii), we observe that by (2.14) we have (again using (2.12)) 



(3.2) 



Tr{e~ tH - e~ tH °' 

oo 
k=2 



)+p ( "\0)t I V(x)dx 



E ,. 



V{X s )ds 



dx 



< 



Pi 



(awv-^iinoc 



'WE 



k=2 



\V{X s )\ds ) dx 

o 



Pi 



C«W,IT,II V-^Woc 



< C7^ ot) (0)t a ||V^||i||^|| oo e 



fc=2 



A:! 



fliv'll 



for some absolute constant C. This proves (ii) and concludes the proof of the theorem. 



□ 



Remark 3.1. We observe that the above holds for any operator which arises from subordination of a 
Brownian motion. For such operators the transition probabilities p t are symmetric radial and decreas- 
ing, hence pt(x, x) = pt(0). In section 5 below we give examples of other nonlocal operators for which 
Theorem 2.2 also holds. 

Proof of Theorem 2.2. We begin by observing that (exactly as in the proof of inequality (3.2)) we have 



e -f*V(x s )ds _ i + T V(X s )ds -\\f V(X s )d 
Jo * Uo 



(3.3) 



\V(X s )\ds, 



for some constant C. By taking expectation of both sides of (3.3) with respect to E x x and then inte- 
grating with respect to x, we obtain 



/ & [ e -foV(x s )ds_ 1+ f V{Xs )ds- l -\fv{X s )d 
jR d V Jo 1 Uo 



n 2 



dx 



< C(t\\V 
= C(t\\V 



I )V l|y|1 

oo I c 



Et 



\V(X s )\ds dx 



) 2 e^ v ^t\\V\\i, 



where we use again (2.12). 

Returning to the definition of the trace differences in (2. 10) we see that this leads to 



1 



pJ q) (o) 



-{Tr{e 



-tH 



(3.4) 



2„t V 



)) + t V(x)dx 



't\\V\U. 



1 







2 \ 


/ E x,x 


[ V(X s )ds 


dx 




Jo 





It remains to estimate the term E x x ([-] 2 ). Since V is uniformly Holder with exponent 7 and constant 
M, we have 



(3.5) 



\V(X s + x)-V(x)\<M\X s \i. 
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Hence, 













E l 


[ V{X s )ds 


2 -t 2 V\x) 




E l 

^x,x 




Jo 









El 



V{X s )ds 



V{X s )ds 



V(x)ds 



V{x)ds 



(3.6) 



By employing (3.5), we obtain 



#o,o ( 



(V(X S + x) - V{x))ds 



I V(X s + x) + V(x))ds 
Jo 









E l ( 

x,x 1 


/ V{X s )ds 


j -t 2 V 2 (x) 




Jo 













< ME% fi ( 


J \X s \^ds 




[I 









(\V(X s + x)\ + \V(x)\)ds 



Integrating both sides of this inequality with respect to x and using Fubini's theorem, we have that the 
value of the second integral becomes 2t||V||i. Thus we arrive at 



(3.7) 



El 



V(X s )ds 
< 2tM||V||i.E§ i0 (J \X s pds\ . 



dx - t 2 



f \v( 



x)\ 2 dx 



Now, it remains to estimate the expectation on the right side of (3.7). As in (2.1 1) we have 



(3.8) E% )0 (J \X B \ids) = J El (\X s y)ds 

r* r pi a) (o,y)pj-i(^o) 

Jo Jri p { t a) (0,0) 
To estimate the right hand side we recall the following "5P-inequality" 



y^dyds. 



(3.9) 



ph t) {x,z)pf 1 {z,y) 



<C a4 (pi a \x,z)+p ( t a) (z,y)), 



valid for all x,y, z £ s, t > and < a < 2, where C a ^ is a constant depending only on a 
and d. This inequality is proved by Bogdan and Jakubowski in [9, p. 182]. It is derived form the the 
"3P-inequality" 



(3.10) 



pi a \x,z)A P r(z,y)<C a4 p^ t (x,y). 



» 



(a) 



(a) 

Inequality (3.10) itself is proved by a simple computation replacing p s with the quantity from (2.7). 
From (3.10) and the fact that for a, b > we have ab = (a A b) (a V b) and (a V b) < (a + &)), (3.9) 
follows immediately. 
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^0,0 



Returning to (3.8), we have, 

" \X a \tds) < C a4 I I ( P i a \0,y)+ P t ] s (y,0))\yrdyds 
J Jo Jm d 

2C a>d I I pi a \0,y)\yrdyds 

Jo JR d 
rt 



2C a>d [ E°(\X s \^ds 
Jo 

2C a . d [ S ^ a E\\X^)ds 
Jo 



2C a4 E°(\X 1 \'V) 



where we used the scaling property of the stable process coming from (2.6). We now recall that 
E°(\Xi\' y ) is finite under our assumption that 7 < a. (This fact follows trivially from (2.7).) Thus 
we see that 

/ f* 

(3.11) 



where the constant C a 7 d depends only on a, 7 and d. We conclude that 











/ E l,x 


f V(X s )ds 


| dx-t 2 f \V(x)\ 2 dx 


< 2tM\\V\\iE% fi 


Jm d \ 


Jo 


J JR d 



\X q Pds 



(3.12) 

Returning to (3.4) we find that 
1 



Tr(e- tH - e~ tHo ' 



< M\\V\\xC a ~ d t^ a+2 . 



))+*/ V(x)dx--t 2 \V(x)\ 2 dx 



(3.13) 



P?(0) " ' 'Ju d ' v ' 2 

< ct 3 ||y||^e t||v||oo ||y||i + M\\v\\ l c an4 v la+2 
<c a , 7 , d ||y|| 1 (|F||^ e *ii 1/ ii-t 3 + ^+ 2 ). 

Rewriting this in the form stated in Theorem 2.2 we arrive at the announced bound 



(Tr(e 



-tH 



))+P?W[ V(x)dx-p { i a \0)t 2 U \V{x) 2 



l dx 



(3.14) <c a , 7id ||y|| lP r ) (o)(|y||Le*ii v 'ii-t 3 + ^ +2 ), 

valid for all t > 0. 



□ 



Remark 3.2. We remark that in the case of Brownian motion, a = 2, one may use the fact (as done in 
[6]) that 



(3.15) 



p?)(o,o) - p t(-) (0 ' y) 



to explicitly compute the right hand side of (3.8). While the inequality (3.9) does not hold for a = 
2, we may estimate the right hand side of (3.8) in the case of Brownian motion without the explicit 
estimate that comes from (3. 15). Indeed, by observing that in law the Brownian bridge started at and 
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conditioned to return to in time t is the same as B s — jBt, where B s is standard Brownian motion, a 
simple computation leads to the same conclusion as (3.1 1) for a = 2. Indeed, 

K,o( [Wds] = [ E°\B s -^B t pds 



o 



S\ 7 



E U [\B s r+[~) \B t p)ds 

C^E 



where C 7i ^ depends only on 7 and d. 

Unfortunately, the simple path construction X s — jXt only leads to the stable bridge when a = 2 
and we are not able to repeat (3.16) in the case < a < 2, hence our use of (3.9). However, it may 
be that the construction by Chaumont in [12] (or perhaps some of the estimates of Fitzsimmons and 
Getoor in [21]) can be used to bypass the estimate (3.9). 

4. Alternative proofs of Theorems 2.1 and 2.2 
Recalling formula (2.14) we have 



k 

)ds I dx, 



(4.1) Tr ( e -*» - e"^) = pM(0) f) U£ f ( f V(X s )ds) 

where as mentioned there the sum is absolutely convergent for all potentials V E L°°(]R ci ) n L 1 (]R d ) 
and all t > 0. We note also that in fact, for every N > 1, 



(4.2) Tr (e~ tH - e- tH °) = (0) ^ ^ V(X s )ds^ dx + 0(t N ^ 



as t J, 0. This follows from the fact that 



(«) E -TT d Ei 4 V{Xs)ds ) M ~ S ^II^H^II^IIi 

II r II 00 1 

(iV + 1) 

f\{X s )d^j =k\J k 



< 1 II ^ Moo II v 111 



Next observe that 



where 

J k = I V{X Sl )d Sl V(X S2 )ds 2 - ■ - V(X Sk )ds k . 

J 0<s l <s 2 <-<s k <t 

Recalling that the finite dimensional distributions of the stable bridge are given by 

1 k 

p lxi x si G dy u X S2 E dy 2 ,...,X Sk E (%J = 

where yo = yfc+i = x, sq = 0, s^+i = t, we can write 
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oo „ 

(4.4) Tr (e- tH - e' tH °) = V(-l) fc / / fc («i, s 2 , . . . , s fc )d Sl d S2 • • • 

fc=1 Jo<si<s 2 <---<Sfe<t 



where 

h(si,s 2 ,. ■ ■ ,s k ) = / Tp^L-s.Cz/j+i -yj)T\V(yj)dyidy 2 .---dy k dx. 

V d ) fc+1 J=0 £=i 

Here (as several times before) we used the fact that pf 1 (x,x) = p[ (0) . By Fubini's theorem we may 
integrate with respect to x first and use the fact that 

/ #(2/1 - x)r {a) ~^-~ (a) 

Jw. d 

to arrive at 



p^(yi - x)p { p Sk (y k - x)dx = pf_! {sk _ Si) {y k - yi) 

oo . 

(4.5) Tr (e~ tH - e"^ ) = ]T(-l) fc / h(si, s 2 ,..., s k )d Sl ds 2 ■ ■ ■ ds k , 

k= l J0<s 1 <S2<---<s k <t 

where 

i k ( Sl ,s 2 ,...,s k ) = j^ k (pt\ Sk . Sl) (yk-yi)v(yi))- 

f k-l \ 

s=i ) 

We can then also re-write the right hand side of (4.2) as 



(4.6) ^ Q) (0) (y^K^- I I k {s 1 ,s 2 ,...,s k )ds 1 ds 2 ---ds k + 0(t I 

\,. ,7)1 '(Q) J0<si<so<~<su<t 



Kk=l P t (0) ■/0<Sl<S 2 <-<8fc<t 

From these formulas we can also compute coefficients in the asymptotic expansion of the trace as 
1 1 fairly directly and in particular give alternative proofs of Theorems 2.1 and 2.2. 
If k = 1, we see that 

pt\s k -s l) (y k -yi)=p ( ? ) (°) 

and 

h(s 1 )=p\ a \0) [ V(y)dy, 

JR d 

which immediately leads to (2.18) in Theorem 2.1. 
If k = 2, we write 

(4.7) h( Sl ,s 2 



+ 




-L 


' f & 

l d JR d 


? (aa _ ai )(ite-yi)p| 


,(<*) 

t-(s 2 - 


- Sl )(^- 


i/i)P(2_ ai) (jte-: 


,(<*) 
t-(s 2 - 






,(<*) 
t-(s 2 - 




2/i)P(2_ Sl) (y2-: 



(s 2 -^l) 1 
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Integrating first with respect to y\ we conclude that 

p£L-,Av2 - yi>S_ sl) (^ - yi)v 2 (y 2 )d m d m = P { t a \o) [ \v( y )\ 2 d y . 

On the other hand, | V{y\) — V(y 2 )\ < M|yi — 2/2 1 7 an hence 



(4.8) 



< M 
= M 
= M 



p£(„_ 81 )(«2 - fi)P(2_ sl )(«2 - Vi) (V(yi) - V{y 2 )) V(y 2 )d yi dy, 

P t _ (S2 _ sl )(y2 - yi)pj„_ ai )(ite - yi)|yi - y2| 7 |^(y2)|dyidy 2 



, ^t— (S2— Si) ^ ■'■*'(. 

From inequality (3.9) we again obtain 



£i 2 _ si) (y)p[l si) (y)|yr^) (jf , |v(y)l<fo) • 



Thus 



where 



(4.9) 



^Vs>)<U)(^r^ ^ c M p^(o) ^ (^ 2 _ si) ( y ) +pg_ ai) (y)) \y\' ( dy- 

I h(si,s 2 )ds 1 ds 2 =p < t a \0)t 2 [ \V( y )\ 2 d y + R t 



\Rt\ < c a4 M\\v\\ lP ^(o) f (£;°|x t _ (S2 _ sl) | 7 + £; |x (S2 _ sl) | 7 )d s 

J0<s 1 <s 2 <t 

< C a ^ d M\\V\\ lP { t a) (0)t^ a + 2 . 



\ds 2 



As in (3.13), this gives another proof of (2.20) in Theorem 2.2. 

5. Extension to other non-local operators 

Taking < j3 < a < 2 and a > 0, we consider the process = X t + alt, where and Y t 
are independent a-stable and /3-stable processes, respectively. This process is called the independent 
sum of the symmetric a-stable process X and the symmetric /3-stable process Y with weight a. The 
infinitesimal generator of Z a is A"/ 2 + a^A^/ 2 . This again is a non-local operator. Acting on functions 



\f(x + y)-f(x)]dy, 



f E Cg°(]Er) we have 



A a/2 + a , A ,/2\ f{x) = A H f /I + 



y 



where Ad, a is as in (2.8). 

Properties of the heat kernel (transition probabilities) for this operator have been studied by sev- 
eral people in recent years. See for example Chen and Kumagai [19], Chen, Kim and Song [17] or 
Jakubowski and Szczypkowski [25], and references given there. If we denote the heat kernel of this 
operator by p%(x), we have that 

f -ix.£ e -t(|f|«"+o< J |f|'') d £ 



(5.1) 
(5.2) 



Pt(x) 



(2ir) d 



1 



(47Ts) d / 2 



e 4s r]f(s)ds, 
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where rjf(s) is be the density function of the sum of the a/2-stable subordinator and a 2 -times the /3/2- 
stable subordinator. Again, this density is radial, symmetric, and decreasing in \x\. It is proved in [19] 
and [ ] that there is a constant C a $^ depending only on the parameters indicated such that for all 

x G R d and t > 0, 

(5-3) C-^J t a (x) < pUx) < C asM f t a (x) 

where 



+ 



\ x \d+a \rj,\d+ft 



We note that for a = this is just the estimate in (2.7). For any 7 > with 0<7</3<awe have 
for any t > 0, 

I E°{\Z a s [<)ds < cJ I E°{\X s yi)ds + a? [ E°{\Y s \^)ds 
Jo \Jo Jo 

(5.4) = CJE\\X^) - ■ - +a^E (\Y 1 y)- 



/ t7 /Q+l J7//S+1 
C a ,a,/3,d — 7 ' 7" + 



7/Q + 1 + I) ' 

It also follows from the "3P-inequality" (24) in [25] (using again the fact that for any a, b > we have 
ab = (aA b){a V b) and (a V 6) < (a + £>)) that we have the following "5P-inequality" for this process: 

f c c\ Ps\ x > z )Pt\ z iV) ^ n (at \ , at w 

(5-5) — < C ajaAd (p B (x, z) + Pt (z, y)) , 

valid for all x,y, z G M. d and all s, t > 0. 

With (5.4) and (5.5), both proofs above for Theorem 2.2 can be repeated to obtain the same result 
for the operator A a / 2 + a 13 A' 3 / 2 , replacing A a / 2 . More precisely we have 

Theorem 5.1. Let a > 0, < /3 < a < 2 an d let Hg = A Q / 2 + a? A' 3 / 2 . Suppose V G L°°(M. d ) n 
L 1 (M. d ) and that it is also uniformly Holder continuous of order 7, with < 7 < (3 A 1, whenever 
< /3 <l,and with < 7 < 1, whenever 1 < /? < 2. Let H a = A Q / 2 + aPkPI 2 + V. Then for all 
t > 0, 

-iff" „-tH%\\ , a/ n \j. / t//„\j„ „a/rV\ ^ j.2 



(Tr(e-* H - e -^o)) +p -(o)t / V (x)dx - p a t (Q)^ 2 \V(x)\ 2 dx 

JR d 1 JR d 

(5-6) < C a>aAjtd \\V\\ lP ?(0) (\V\\ 2 co e t ^t 3 + t'y/ a+2 + f^+ 2 ) , 

where the constant C a ^ a ^ n ^ depends only on a, a, (3, 7 awe? d particular, 

(5.7) Tr{e~ tHa - e~ tH S) = Pt a (0) (-4 / F(x)dr + \t 2 [ \V{x)\ 2 dx + 0(t^ a+2 )) , 

V JR d ^ Jjgd J 

as t 

We next consider another stochastic process associated with a non-local operator that has been of 
considerable interest to many people in recent years, the so called relativistic Brownian motion and its 
more general versions, the a-stable relativistic process. This is again a Levy process denote by X™ 
with characteristic function 
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(5.8) e -t^e +m ^ 2 -m) = E{S .x r) = f e *-ypt' a) (v)dy, 

for any m > and < a < 2. The case a = 1 is the relativistic Brownian motion and m = yields 
back the stable processes of order a. Henceforth we assume that m > 0. As in the case of stable 
processes, X™ is a subordination of Brownian motion and in fact 

(5.9) P [ r a \*)= = f e-<e-i^ 2+m2/ ^ /2 - m )d^ 

: e 4s r] t [s) as, 



( 47rs )d/2 

where rff 1 ' 01 ^) is the density of the subordinator with Bernstein function $(A) = (A + m 2 / a ) a l 2 — m. 

As before we see that p[ m ' a \x) is a radial, symmetric, and decreasing in \x\. By a simple change of 
variables it is easy to see that the transition densities have the scaling property 

(5. 10) p[ m ' a) (x) = m d l a £f ) (m^x) . 

These properties can all be found in Ryznar [3 0] where it is also proved that 

(5-11) pt' a) (x) = e mt r - \ f!/9 e^ 1 e- m2/as r ] f 2 ( S )d S , 



lo (47Ts) d / 2 

where rj^ 2 (s) is the density for the a/2-stable subordinator. By scaling 

11 f 2 {s)=t-^\st-^). 
Hence changing variables in (5.1 1) leads to 



(5.12) lime- mt t d / a pl m ' a \0) = / rft /2 {s) ds = p?(0) 



40 n v ' J (47rs) d / 2 11 " J J (2vr) d a ' 

where the last equality follows from (2.5). 

The infinitesimal generator of X™ is given by m — (m 2 ' a — Af 2 . The case a = 1 gives the 
generator m — V— A + m? which is the free relativistic Hamiltonian. For more on these notions we 
refer the reader to Carmona, Masters and Simon [11]. 

Estimates for the global transition probabilities pj m,Q ^ (x) and their Dirichlet counterparts for various 
domains have been studied by many authors. We refer the reader to [13], [14], [15], [16], [18], [19], 
[29], [30]. Here we recall the following small time estimate from Chen, Kim and Song [ , Theorem 
2.1]. There exists a constant C a%m% d depending on a, m and d such that for all iel d and all t £ (0, 1], 

(5-13) C-'-* A < p <-> (I) < (^f* A 



where 



* (r) = 2 -(^) r fd + °\ 1 a ^-l e -'/4 e -r a /* da 



which is a decreasing function of r with ^(O) = 1, \&(r) < 1 and with 

cr l e -r r (d+a-l)/2 < ^ (r) < Cie - v ( d+a -l)/ 2) 



14 



RODRIGO BANUELOS AND SELMA YILDIRIM YOLCU 



for all r > 1. (See also Chen and Song,[18, p. 277].) Setting 

pt(x) = A 



I 

it follows trivially (exactly as in Bogdan and Jakubowski [9, p. 182]) that 
(5.14) p s (x) Apt-six) < Cp t (0), 

for all iel d and < s < t < 1. Therefore, 

(5-15) p^ a) (0,x)Apt' s a \x,0) < C a ^ d p { r a) (0,0) 

for all x G K d and < s < t < 1. This as before leads to the "5P-inequality" 

(5-16) ^^^V^ < <W fe m ' a) (x,0) +^(0,x)) , 

valid for all x G M. d , < <s < t < 1. (Recall that pi m ' Q ^(x, = pi m ' Q ^(x — y).) Even though this 
estimate is not as general as the previous "5P-inequalities" above, it suffices for our needs. 

Next, we need to compute bounds for the 7 moments of X™. More precisely we want to understand 
the behavior of f Q E°\X™\" t dr, as t J, 0. Let us denote the standard Brownian motion in M d by B s and 
recall that E°\B 2s \' r = T^^s^^E^B^ and 2t/ 2 £°| j Bi|t = C ljd , where the latter is an explicit and 
computable constant whose exact value is not relevant for our calculations here. Using this, (5.10) and 
(5.11) we have 

t , ft rmt 

E°\X^ds = / E°\Xi r \^dr = m ( ^~ 1} / E°\X l T [<dr 

Jo Jo 

(5.17) = m ( ^ _1) / / \x\ n/ p^(x)dxdr 

(±L-1) 



rmt / roo \ 

J e r I J E \B 2s \^e- s 7]^ 2 (s)dsjdr 

rmt / roo \ 

C^ d m^-V J e r (J s^/ 2 e- s ^l 2 {s)dsjdr, 
C^m^-V J™' e r QH s^ 2 e- s r- 2 / a ^ /2 (r- 2 / a s)ds^ dr 
C^m^J^r^ (jH s<l 2 e~ sr %? /2 ( S )<fe) dr, 



where we used the scaling property of the stable subordinator for the second to the last line and a 
change variables for the last line. But 

roo roo 

/ s^ 2 e- sr2/a ^ /2 (s)ds < / s^X /2 (s)ds = C J E°\X 1 p < 00, 
Jo Jo 
whenever 7 < a. (Here, X\ is the stable process of order a at time 1.) Thus, 

f E°\X?rdr < C a ^ nhd t^ a+1 e mt 
Jo 

and this is valid for all t > 0. From this and the inequality (5.16), we obtain 
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£0.0 (flX^ds) < C a , m4 f f (p^ a) (y,0)+ptf\0,y))\yPdyds 
\Jo J Jo Jm d 

(5-18) = 2C a , m . d I I p^ a \o,y)\y^dyds 

Jo J~R d 

= 2C aimA [ E°(\X™F)ds 
Jo 

_ ^0,7,771,(1 b J 

for all t G (0,1]. 

With (5.18), the proof of Theorem 2.2 can be repeated to obtain a similar result for the operator 
associated with the relativistic stable processes. More precisely we have 

Theorem 5.2. Let H™ = m - {m 2 l a - A) a/2 . Suppose V G L°°(R d ) n L 1 (R d ) and that it is 
also uniformly Holder continuous of order 7, with < 7 < a A 1, whenever < a < 1, arc J w;7/? 

< 7 < 1, whenever 1 < a < 2. Let H m =m - (m 2 / a - A) a/2 + V. Then, for all t G (0, 1], 
(Tr(e- tHm -e- tH ™))+pU0)t [ V(x)dx - p?(0)^ / |F(x)| 2 dx 

(5.19) < C^.llFHi^^O) (|yf ooe *^ll-t 3 + ^+ 2 ) , 

where the constant C arurn ^d depends only on a, 7, m and d. 
In particular, 

(5.20) Tr(e- tHm -e- tH ™)=p ( r> a) (0)(-t [ V{x)dx + \t 2 [ \V(x)\ 2 dx + 0(t^ a+2 ) ), 

\ jR d 2 y K d / 

i 4, 0. 
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